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A route to room temperature ferromagnetic ultrathin SrRu0 3 films 
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Experimental efforts to stabilize ferromagnetism in ultrathin films of transition metal oxides have 
so far failed, despite expectations based on density functional theory (DFT) and DFT+U. Here, 
we investigate one of the most promising materials, SrRuOs, and include correlation effects beyond 
DFT by means of dynamical mean field theory. In agreement with experiment we find an intrin¬ 
sic thickness limitation for metallic ferromagnetism in SrRu 03 thin films. Indeed, we demonstrate 
that the realization of ultrathin ferromagnetic films is out of reach of standard thin-film techniques. 
Proposing charge carrier doping as a new route to manipulate thin films, we predict room temper¬ 
ature ferromagnetism in electron-doped SrRuOs ultra thin films. 


PACS numbers: 31.15.A-, 31.15.V-, 73.50.-h, 73.61.-r 

Thin films and heterostructures of the 4 d perovskite 
SrRuOs (SRO) are intensively studied and used, in par¬ 
ticular, as gate electrodes for novel oxide-based electronic 
devices [1, 2]. The reason for this is that SRO is a conduc¬ 
tor with good thermal properties [3] and high chemical 
stability, allowing for epitaxial growth on various sub¬ 
strates, as well as the combining with other perovskite- 
based materials to form complex heterostructures [4, 5]. 
In the bulk, SRO is a ferromagnetic (FM) metal with 
a, for a 4 d oxide, remarkably high Curie temperature, 
Tc = 160 K, and an experimental magnetic moment 
in the range of 0.8 to 1.6 t l B [6-9]. SRO further at¬ 
tracted fundamental research interests regarding, among 
others, magnetic monopoles [10], non-Fermi Liquid [11], 
spin freezing [12], and the debate of itinerant [13] versus 
localized magnetism [14]. 

However, the FM moment and Curie temperature get 
dramatically suppressed below a sample thickness of 4 
unit cells [15-18], and eventually single unit cell SRO 
films turn antiferromagnetic (AF) and insulating [17, 18]. 
This led to the pertinent question whether there is a 
fundamental thickness limit for ferromagnetism [18], and 
concerted efforts to stabilize ferromagnetism in ultrathin 
SRO films by compressive and tensile strain or capping 
layers [19, 20]. However, hitherto ferromagnetism in sin¬ 
gle unit cell films remains unattainable for SRO or any 
other oxide material, even in a heterostructured setup. 

On the theoretical side, previous attempts to under¬ 
stand the electronic structure and the transition to an 
AF insulator resorted to density functional theory (DFT) 
and the static mean-field DFT+U approach. The for¬ 
mer failed to reproduce the transition [21], while the lat¬ 
ter found a transition to an AF insulating state below 
four layers when assuming an artificial Ru 02 terminated 
surface [22], while experimentally samples are found to 
have a SrO termination [5]. DFT+U further predicted a 
spin-polarized highly confined half-metallic state for an 
SRO mono-layer when either sandwiched with SrTiOs 
(STO) [19] or grown on a strained STO substrate [20]. 
However, such a state could not be confirmed in exper¬ 
iment [23]. The apparent discrepancy between experi¬ 


ments and results from standard band-structure methods 
calls for a more sophisticated treatment of electronic cor¬ 
relation effects. Indeed already in the bulk, SRO displays 
signatures of electronic correlations, such as many-body 
satellites in photoemission or the violation of the Ioffe- 
Regel limit in the resistivity [8, 24]. Hence, SRO is to be 
considered an -at least- moderately correlated system. 
Note that a dimensional reduction/geometric constraints 
in thin films can be expected to further enhance elec¬ 
tronic correlations. 

For a better and unbiased treatment of these corre¬ 
lations effects in various SRO films and heterostructure 
setups, we employ realistic DFT + dynamical mean-field 
theory (DMFT) [25-29] calculations. Our main findings 
are: (1) Both the SRO mono-layer and bi-layer are AF 
insulators. (2) We demonstrate that standard thin film 
manipulation techniques such as strain and surface cap¬ 
ping can neither restore ferromagnetism nor metallicity 
to a SRO mono-layer; interestingly, we find that surface 
capping pushes the AF insulator towards a paramagnetic 
(PM) insulator. (3) With new insight regarding the mi¬ 
croscopic origin for the transition, we identify carrier dop¬ 
ing as the best option to generate FM properties that 
are on a par with those of the bulk. We find the FM 
moments of doped SRO films to be stable even at room 
temperature, heralding a great potential for technological 
applications. 

Method. We use the experimental orthorhombic crys¬ 
tal structure of SRO [30] for the various setups of bulk, 
films and heterostructures. In the films and heterostruc¬ 
ture both the internal positions and lattice constant are 
relaxed; the in-plane lattice constants in films are fixed 
to the experimental ones of STO. Fig. 1 exemplary shows 
the SRO mono-layer grown on 4 layers of STO sub¬ 
strate. The atomic relaxations are carried out with the 
VASP program package [31, 32] using the PBE functional 
[33]. For the optimized atomic positions, we subsequently 
perform WIEN2K [34] electronic structure calculations 
with the mBJ exchange [35] and PBE correlation func¬ 
tional [36], and a Wannier function projection onto max¬ 
imally localized [37] t 2 g Wannier orbitals [38] using the 
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Wien2Wannier program package [39]. This t^g Hamil¬ 
tonian is supplemented by a local Kanamori interaction 
and solved within DMFT using Wien2Dynamics [40] , em¬ 
ploying a hybridization expansion [41] continuous-time 
quantum Monte Carlo (CTQMC) algorithm. For the 
Coulomb interaction strengths, we adopt a Hund’s ex¬ 
change (J = 0.3eV), intra- (U = 3.0eV) and inter¬ 
orbital Coulomb repulsion (U' = 2.4eV). These values 
are chosen not only because they are in between the con¬ 
strained random phase approximation (cRPA) values for 
(i) free standing cubic SRO mono-layer (0.3eV, 3.5eV, 
and 2.9eV) and (ii) orthorhombic bulk and (0.3eV, 2.3eV, 
and 1.7eV), but also because they reproduce the FM 
metallic state for orthorhombic bulk SRO. 

Bulk SrRuOs. The moderately correlated electronic 
structure of bulk SRO was successfully captured in both 
many-body perturbation theory [42] and realistic DMFT 
calculations [43, 44]. Also DFT calculations correctly 
predict that SRO is an itinerant ferromagnet with a mo¬ 
ment ranging from 1.5 to 1.6/ib [8, 45]; similar moments 
have also been obtained within DFT+DMFT [44, 46]. 
Using DFT+DMFT, we indeed find orthorhombic bulk 
SRO to be a FM metal with orbital occupations of 0.867 
(0.466) for the majority (minority) spin of all three t^g or¬ 
bitals at T = 100K. This corresponds to a FM magnetic 
moment of ~ 1.2 hb- A GdFeOs-type distortion in which 
the corner-sharing octahedral tilt around the y -axis and 
rotate around the 2 -axis lifts, in principle, the t^g degen¬ 
eracy. The effect on the crystal field and orbital occu¬ 
pations of orthorhombic SRO is however minute. Our 
DFT+DMFT finds SRO to be a PM metal above T c ~ 
150K which is close to the experimental Curie tempera¬ 
ture of 160K (cf. Fig. 3 below). 

Thin films. We now consider SRO grown on STO, and 
study the evolution of the electronic structure when re¬ 
ducing the number of SRO layers. We find that FM is 
suppressed: the mono- and bi-layer SRO on STO are AF 
insulators, in congruence with experiments [18, 23] that 
show a dramatic drop in the FM moment and an insulat¬ 
ing behavior for 4 layers or lower [15-18, 23]. Indicative of 
an itinerant origin of ferromagnetism, the critical thick¬ 
ness of the magnetic and electronic transition coincide. 
Fig. 2 (a) shows the spectral function of the mono-layer. 
The system is gapped by ~1.0eV and displays a large 
orbital polarization: the xy orbital is fully filled, and the 
xy and xz orbitals are half-filled and fully spin-polarized, 
resulting in an AF moment of ^ 2 /is- This finding is 
supported by recent exchange bias measurements [17]. 

We note that for the particular case of the SRO mono- 
layer, also LDA+U [20] can seemingly give a qualita¬ 
tive correct picture, as the system is orbitally and spin- 
polarized. However, the underlying physics is very differ¬ 
ent: When heating the mono-layer above its Neel tem¬ 
perature within DMFT, the system remains insulating 
at non-integer filling [0.88 (0.88), 0.56 (0.56), and 0.56 
(0.56) for the spin up (down) xy , yz and xz orbitals at 



FIG. 1. Right: Structure of a SRO mono-layer grown on 4 
layers of STO. Upper left: Top view of the same structure. 
The indicated y/2 x y/2 supercell was adopted for allowing 
AF-ordering in each Ru02 layer. Lower left: atomic labels 
and coordinate system (Figures drawn with the Vesta code 
[47]). 



FIG. 2. DFT+DMFT spectral functions of (a) the AF SRO 
mono-layer on a STO substrate and (b) the FM doped (4.3 
electrons/Ru) superlattice (STO)s:(SRO)i at 150K. Insets: 
electronic occupations. 


1000K]. This complex Mott physics [48] reveals that the 
AF insulating phase is beyond a simple Slater descrip¬ 
tion, and thus not describable by LDA+U. 

Physical origin of transition. Let us now investigate 
the microscopic origin of the FM-metal to AF-insulator 
transition. Whereas the crystal field splitting is minute 
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TABLE I. Magnetic and conductive states of bulk SRO, 
SRO mono-layer (also under 0.5% compression and 0.5% 
tension), SRO bi-layers and (STO) 5 :(SRO)i superlattice 
(capped mono-layer), as obtained from spin-polarized DFT 
and DFT+DMFT calculations in comparison with experi¬ 
ment. FM-M: ferromagnetic metal, AF-I: antiferromagnetic 
insulator; NM-I: non-ferromagnetic insulator (the magnetic 
nature of the experimental non-ferromagnetic state has not 
been fully determined yet; the exchange bias behavior hints 
at antiferromagnetism [17]). 


System 

DFT 

DMFT 

EXP 

bulk 

FM-M 

FM-M 

FM-M [ 6 ] 

mono-layer 100.5% 

FM-M 

AF-I 

- 

mono-layer 100 % 

FM-M 

AF-I 

NM-I [17] 

mono-layer 99.5% 

FM-M 

AF-I 

- 

bi-layer 

FM-M 

AF-I 

NM-I [15, 16] 

superlattice 

FM-M 

AF-I 

NM-I [23] 

hole doped superlattice 

FM-M 

FM-M 

- 

electrons doped superlattice 

FM-M 

FM-M 

- 


for the bulk, in the SRO mono-layer the x?/-orbital is 
energetically lower than the yz and xz orbitals, because 
the (cubic) crystal symmetry is strongly broken at the 
surface. This is already the case for the DFT Wannier 
Hamiltonian, but correlation effects boost the crystal 
field splitting [48-50] of the SRO mono-layer (see Sup¬ 
plemental materials Table I). Therewith, the xy-oibitol 
become essentially fully occupied, and the two remain¬ 
ing electrons occupy the yz and xz orbitals: The single¬ 
layer SRO is an effective half-filled two-band system, 
favourable to AF order. 

Tuning the properties of the SRO mono-layer. The 
prime motivation for SRO-based thin films are the advan¬ 
tageous properties of the FM metallic bulk. However, the 
desired features such as the magnetic moment strongly 
decrease for thinner films and eventually ultrathin films 
become non-FM, in agreement with our calculations. A 
natural question is whether the bulk properties can be 
restored, at least partially, by tuning the geometry of the 
films. 

First we discuss the influence of straining/tensioning 
the mono-layer. This can be realized experimen¬ 
tally by choosing an appropriate substrate. Indeed, 
previous DFT+U calculations [20] predicted a strain- 
induced FM half-metallic state for the SRO mono-layer. 
DFT+DMFT, however, does not show any tendency to¬ 
wards a FM half-metallic state at least for realistic U 
values, see Table II in the Supplemental Materials. Also, 
the effective crystal-field splitting A e ff, shown in Table 
I of Supplemental Materials, can only be tuned slightly 
through straining/tensioning. 

Another way to influence the crystal-field splitting is 
through the deployment of capping layers. Here, we 
study the effect of capping the SRO mono-layer with 
additional layers of STO. Specifically, we consider a 
(STO) 5 :(SRO)i superlattice [19] consisting of 5 layers of 
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FIG. 3. (a) Magnetic moment of SRO mono-layer (ML) and 
(STO) 5 : (SRO) i superlattice (SL) vs. doping at 150K and 
300K (positive [negative] moments denote FM [AF] ordering), 
(b) Magnetic moments vs. temperature for orthorhombic bulk 
SRO, mono-layer, undoped and doped (STO) 5 :(SRO) 1 super¬ 
lattice. 


STO alternating with a mono-layer of SRO. This restores, 
at least partially, the hopping amplitudes in the out-of- 
plane direction. Compared to the SRO mono-layer, the 
DFT crystal field splitting, shown in Supplemental mate¬ 
rials Table I, is now much smaller (-0.05eV) approaching 
the negligible value of the bulk. As a result, the t^g 
orbital occupations are more balanced. However, this 
causes only a slight reduction of the AF magnetic mo¬ 
ment (1.92 fiB at 150K and 1.48/i# at 300K, cf. Fig. 3 
(b)) with respect to the un-capped mono-layer. Our find¬ 
ing of a non-FM insulating state with a gap of ~1.0eV 
for the capped mono-layer is consistent with experiments 
[23, 51], where it was concluded that SRO capped by STO 
leads to a insulator without a net moment. We note that 
previous DFT+U calculations [19] instead predicted a 
FM half-metal, at variance with experiment. 

The above calculations reveal that the non-FM insu¬ 
lating state is a robust feature of the SRO mono-layer. 
Only for an unrealistically small U-to-bandwidth ratio a 
FM phase can be stabilized (see Supplemental materi¬ 
als). In consequence, none of the standard manipulation 
strategies available to the production of thin films can 
tune this ratio sufficiently to induce ferromagnetism to 
the SRO mono-layer. This stability explains why all ex¬ 
perimental efforts to create ultrathin ferromagnetic films 
have so far been unsuccessful. 
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Doping. Here we propose an alternative route to 
achieve ultrathin FM films: doping. This strategy may 
seem counter-intuitive at first glance, as carrier doping 
causes a deterioration of the desired properties in the 
bulk. As we shall see, the situation for ultra thin films 
is different: Using the virtual crystal approximation to 
simulate carrier doping within DFT-fDMFT, we obtain 
for the SRO mono-layer and the (STO) 5 :(SRO)i super¬ 
lattice the magnetic moments shown in Fig. 3. For the 
SRO mono-layer a significant doping corresponding to 3.5 
and 4.5 electrons/site is needed to turn the AF state into 
a PM at both low (150K) and high temperature (300K), 
see Fig. 3 (a). The AF magnetic moment is essentially 
symmetric around the filling with four electrons. The 
reason for this is the previously mentioned crystal field 
effect that results in an almost fully occupied xy orbital 
and a half-filled (particle-hole symmetric) yz/xz orbital 
doublet. We note that in all cases, the doping away from 
four electrons induces a metallic state. 

Let us now turn to the more important case, the 
(STO) 5 :(SRO)i superlattice: At low temperatures, e.g. 
at 150K as shown in Fig. 3 (a), both hole and electron 
doping can induce strong FM states. In the case of elec¬ 
tron doping (filling > 4) the FM state is accompanied 
by an alternating orbital ordering of the xz/yz minority 
spin. The spectral function corresponding to 4.3 elec¬ 
trons/site is shown in Fig. 2 (b). There, one Ru site 
has the orbital occupations for up (down) spin: xy 1.00 
(0.34), yz 0.97 (0.83), xz 0.97 (0.19); while for the second 
Ru site: xy 1.00 (0.34), yz 0.97 (0.19), xz 0.97 (0.83). 
For hole doping, e.g. at 3.7 electrons, our DMFT re¬ 
sults indicate that the system is a FM half-metal with a 
moment of 2.20 fis /Ru at 150K, see Fig. 3 (b) and the 
Supplemental Materials for the corresponding spectral 
functions. The half-metallic behavior makes this setup 
a prospective candidate for spintronics applications. To 
put this finding into perspective, we recall that bulk SRO 
has an FM moment of 2/i#/Ru or less, and an experi¬ 
mental (theoretical) Curie temperature of “only” 160K 
(150K). One might thus wonder whether the ferromag¬ 
netism of the doped supercell is actually superior to the 
hailed properties of stoichiometric bulk SRO, which we 
were striving to restore. To investigate this, we perform 
calculations as a function of temperature, see Fig. 3 (b). 
We find that for both hole and electron doping, mag¬ 
netic moments and Curie temperatures of the supercell 
are remarkably higher than for bulk SRO. In particular, 
for 4.2 electrons/Ru site, a sizable magnetic moments 
survives up to room temperature 300K, see Fig. 3 (b). 
The magnetization curve has a similar shape as for the 
bulk [46], despite the much higher T c and the orbital¬ 
ordering. In the Supplemental Material we also go be¬ 
yond the virtual crystal approximation and show that a 
(STO)5:(Lao.25Sro.75Ru03)i superlattice is indeed FM. 
Our findings pave the road for realizing FM oxide de¬ 
vices that can be operated at room temperature. 


Conclusion. Including many-body effects by means 
of DFT+DMFT, we show that the SRO mono-layer is 
an AF Mott insulator owing to a correlation enhanced 
crystal-field splitting at the surface/inter face. While the 
bare (one-particle) crystal-field splitting can be tuned to 
almost zero by STO capping layers, electronic correla¬ 
tions are still strong enough to boost the orbital separa¬ 
tion so that also the capped SRO layer remains an AF 
insulator. A FM metallic state is only realized for an 
interaction-to-bandwidth ratio that cannot be realized 
by experiment. This explains why no ultrathin FM films 
could be stabilized in experiment to date. Given the ro¬ 
bustness of the AF state of SRO mono-layer setups to 
standard thin film manipulation techniques, we propose 
an alternative route to realize a FM state: Our study 
suggests that carrier doping drives ultrathin SRO films 
capped with STO into a strong FM state, whose ordered 
moment and Curie temperature even exceed the values 
realized in stoichiometric bulk SRO. To achieve the long¬ 
standing quest for a FM ultrathin film in practice, we 
consider inducing oxygen or Sr vacancies [52] or doping 
potassium into STO:SRO superlattices [23] as the most 
promising means. 

Our study also opens a new, general, perspective: Pro¬ 
ducing heterostructures based on materials with opti¬ 
mized bulk properties (e.g. stoichiometric SRO) is actu¬ 
ally not always the optimal way for achieving those prop¬ 
erties in a film geometry. Indeed the electronic structure 
of the thin film is so different from the bulk that it can 
be viewed as a completely different material. A manipu¬ 
lation (in our case doping) that decreases the quality of 
the bulk, can in fact enhance the sought-after property 
(FM magnetic moment) for the film setup. This suggests 
in turn that rather inconspicuous bulk materials might 
actually be good candidates for specific functionalities 
when deployed in a film or heterostructure. With this ob¬ 
servation the repertoire of materials to be evaluated for 
oxide-electronics applications is significantly enlarged. 
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